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As a test of quenched lattice twisted mass QCD, we compute the non-perturbatively 0(a) improved pseudoscalar 
and vector meson masses and the pseudoscalar decay constant down to Mps/My — 0.467(13) at /3 = 6 in large 
volume. We check the absence of exceptional configurations and - by further data at /3 = 6.2 - the size of scaling 
violations. The CPU time cost for reaching a given accuracy is close to that with ordinary Wilson quarks at 
Mps/Mv ~ 0.6 and grows smoothly as Afps/Mv decreases. 



1. INTRODUCTION 



The lattice twisted mass QCD (tmQCD) action 
|l] in infinite space-time volume reads 



Am = 7V - f V*V + mo + iMq75T^ 



(1) 



where V = (V + V*)/2, is a flavour quark dou- 
blet and the matrix acts in the flavour space. 
If niQ = nic, the quark mass is entirely given by 
the multiplicatively renormalized parameter /iq. 

Lattice tmQCD is a regularization of QCD with 
N{ — 2 mass degenerate quarks [Q where an ac- 
tion term that is proportional to /iq and is not 
aligned in flavour chiral space with the Wilson 
term contributes to the soft breaking of the chi- 
ral symmetry. As a consequence, the lowest eigen- 
value of the Hermitean square of the Dirac matrix 
on any gauge background is bounded from below 
- in lattice units - by {afiq)^. The problem of the 
exceptional configurations in quenched and par- 
tially quenched QCD is thus solved as long as the 
quark mass is not vanishingly small. 

*Based on a poster presented by R. Frezzotti at the XIX 
International Symposium on Lattice Field Theory "Lattice 
2001", August 19-24, 2001, Berlin, Germany. 



One can think of the two mass-degenerate 
quark flavours as the u and d quarks in a chi- 
rally twisted basis, with the obvious consequences 
on the physical meaning of the correlation func- 
tions 1^ . Heavier flavours of Wilson quarks can 
still be added e.g. in the usual way, while the 
neglected tiny mass difference between u and d 
quarks can safely be taken into account by means 
of chiral perturbation theory. Moreover, lattice 
tmQCD can also be formulated for a doublet of 
non-degenerate quarks, whilst retaining the pro- 
tection against exceptional on all gauge back- 
grounds Q. 

2. OBSERVABLES AND SIMULATIONS 

In this study we restrict attention to the pseu- 
doscalar and vector meson masses, Mps and My, 
the pseudoscalar decay constant, -Fps, and the 
renormalized quark masses, /^r and mpt, which 
are defined according to eq. (3.3) and eq. (3.20) 
of Ref . . It is convenient for our discussion to 
define polar quark mass "coordinates" via 



tan a = fiB,/ , = (/i^ + : 



(2) 



We work in the quenched approximation con- 
sidering systems of physical size L^T at /3 = 6 and 
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Set 


/3 




L/a, T/L 


^ mcas. 


MrTo 


Mpsro 






Al 


6 


0.135196, 0.0266 


16, 2 


650 


0.2729(15) 


1.711(7) 


0.455(5) 


2.662(40) 


Al' 


6 


0.135196, 0.0266 


16, 3 


650 


0.2729(15) 


1.714(6) 


0.455(6) 


2.656(42) 


A2 


6.2 


0.135814, 0.0180 


24, 2 


535 


0.2558(16) 


1.623(8) 


0.456(5) 


2.557(32) 


BIL 


6 


0.135208, 0.0190 


24, 2 


260 


0.1949(11) 


1.452(6) 


0.432(6) 


2.517(35) 


Bl 


6 


0.135208, 0.0190 


16, 2 


535 


0.1949(11) 


1.455(8) 


0.428(5) 


2.513(47) 


B2 


6.2 


0.135814, 0.0138 


24, 2 


300 


0.1962(12) 


1.420(9) 


0.436(7) 


2.462(41) 


C 


6 


0.135208, 0.0117 


24, 2 


260 


0.1205(7) 


1.160(6) 


0.401(6) 


2.485(59) 



Table 1 



(3 = 6.2, with L such that MpsL > 4.5 to suppress 
finite volume effects. Following closely Ref. ||], 
we adopt Schrodingcr functional (SF) boundary 
conditions for tmQCD and compute the SF cor- 
relators 



/A^(xo), f^\xo), fVixo), n 



rll 



(xo) , kj^ (xo) 



(3) 



as a function of the Euclidean time xq , for L rang- 
ing from 1.5 to 2.2 fm, depending on the quark 
mass values, and T/L = 2 or 3. We hence con- 
struct the renormalized and 0(a) improved SF 
correlators [/A']R(a;o), [/p^]R(a;o), [ky,]n{.xo) and 
[fcx^]R(xo), which are specified in eq. (3.11) of 
Ref. 1^ with a given by eq. (^). 

In the limit of large Xq and large T — xq and up 
to cutoff effects, the correlators [/a'(^o)]r and 
[/p}(xo)]r, which correspond to the insertion of 
the operators (A^)o and (P^)^, see eq. (3.10) of 
Ref. ^ , are expected to be dominated by the low- 
est isotriplet pseudoscalar state. In the same limit 
the correlators [fcv/(2;o)]R and [fc^^(xo)]R, which 
correspond to the insertion of the operators (V-^)}. 
and (?R)feoi (3.10) of Ref. are domi- 

nated by the lowest isotriplet vector state. 

In our simulations we employ the non- 
perturbatively improved version of the action (|^) , 
with Csw taken from Ref. . An overview of our 
simulation parameters, statistics and preliminary 
results is given in Table 1, while a = 7r/2 + 0(a). 

Our most critical simulation (set C) required 
^ 230 GFlops X day by employing a CGNE solver 
for the SSOR-preconditioned Dirac matrix Dtm- 
For this (and similar) choice(s) of parameters the 
BiCG solver often does not converge, while a 
CGNE solver works fine with small fiuctuations 
in the number of iterations. 
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Figure 1. The square root of the relative a priori 
variance of / = /p^(24a) versus the number of 
measurements N for the simulation C (filled cir- 
cles) and a simulation with the same values of /3 
and Mr, but /Xq = (open squares). 

3. RESULTS 

We find that lattice tmQCD allows, as ex- 
pected, to safely work in a region of parameters 
which would be inaccessible with ordinary Wilson 
quarks: see e.g. Fig. |l|. For a given number of 
independent measurements, the statistical errors 
on Mps and My are comparable, up to a factor 
of one to three, to those found e.g. with domain 
wall quarks The CPU time effort, e.g. for 
the data sets Al, Al' and A2, is in line with the 
computational cost for ordinary Wilson quarks. 

Concerning the extraction of hadron masses 
and matrix elements, we recall that the lattice ac- 
tion (H) enjoys reduced parity and isospin symme- 
tries I ij: it just preserves the third isospin genera- 
tor and parity combined with a flavour exchange. 
At finite a the correlators may hence receive con- 
tributions of order a^q from the allowed states 
with different parity and isospin. While deferring 
details to Ref. ||], we show in Fig. || an example 
of effective masses extracted from SF correlators: 
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Figure 2. Pseudoscalar effective masses extracted 
from the correlators [/a']r and [/p?^]r for data set 
C. The circles denote our fit to [/a']r: which ac- 
counts for 0(a/iq) contributions from states with 
twisted parity and isospin at large xq. 



the correlator [/^']R(a;o) receives contributions of 
order a/iq that are peculiar to our SF setup and 
would be absent at a = 0. Analogous effects in 
[/p^]r(2^o) are expected to be very small 

As detailed in Ref. we expect the relations 
among our observables and the renormalized pa- 
rameters ro and Af^ to be 0(a) improved. In par- 
ticular, when working at a = 7r/2-|-0(a), an 0(a) 
improved estimate of _Fps is obtained in terms of 
Mps and the matrix element of e.g. the operator 
2Mr(P;^)-'^ between the charged pion and the vac- 
uum[y|. This approach to -Fps, which is also valid 
in unquenched lattice tmQCD, requires to know 
no other counterterms than Zg, rric and Cgw 

In order to check for the residual scaling viola- 
tions, we produced data at (5 — 6.2 (sets A2 and 
B2), while keeping a, Mr and tq fixed. More pre- 
cisely, the ratio m^/AIji was tunecj^ to the value 
0.000(1) for Mnro = 0.195(1) and -0.015(1) 
for Mftro = 0.273(2). For the simulation C 
we found m^/M^ = 0.083(5). We also reanal- 
ysed the data of Ref. 0, which were produced 
at /3 = 6,6.1,6.2,6.45, by imposing precisely 
the same renormalization conditions as in this 
study of tmQCD. We then performed a contin- 
uum extrapolation, assuming a purely quadratic 
dependence on (a/ro)^ and discarding the data at 
P — 6. However, owing to the use of the one-loop 
estimate of 6a, the resulting estimate of Fpg is 



•^The small mismatch in Mj^ro for the set A2 was cor- 
rected by employing estimates of the dependence of our 
observables on Al^rg. 
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* ALPHA + UKQCD P=6. 
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Figure 3. Mpsro and Fpsro versus A/pro from 
tmQCD and reanalysis of the data of Ref. 0, 
including a continuum extrapolation (c. 1.). 

not fully 0(a) improved. The outcome of this ex- 
ercise is compared with the results from tmQCD 
in Fig. ||. We omit the case of Myro, as cutoff 
effects are hardly visible within statistical errors. 

4. CONCLUSIONS 

Lattice tmQCD solves the problem of excep- 
tional configurations and allows for rather precise 
results with a moderate computational effort and 
tiny scaling violations, at least for the observables 
of this study. The framework is suitable for ex- 
tension beyond the quenched approximation. 

This work is part of the ALPHA Collaboration 
research programme. 
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